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1. Introduction 

It is well known from solid state theory [[U 13, x-ray defraction [3J, or the atomic motion 
in laser fields pi] that the presence of a periodic potential requires a modification of the 
conventional concepts of scattering theory In the context of ultra-cold quantum gases, 
this has led to a tremendous outburst of activities, during the past years. Today, it is possible 
to examine the interplay of many-body physics at the lowest attainable temperature, in the 
presence of designable optical lattices O for bosons EH HI [Hi, fermions lITTllTIl fBlfTill. 
bose-fermi mixtures lfT5l[T6l and to manipulate simultaneously the interaction among particles 
|[I71[l8l[l9ll20ll2ll, as well as collective states ri22il . 

In the present article, we will discuss binary Feshbach resonance scattering in the 
presence of a lattice, as a particular aspect of the aforementioned general theme. When one 
approaches the topic of binary atomic scattering in homogeneous space and in a periodic 
lattice, one needs to highlight the similarities and differences, first. In homogeneous space the 
scattering process connects asymptotical free states, which are plane waves \q) . The accessible 
relative kinetic energy > forms a simple continuum, which has a lower bound, but 
no upper one. Furthermore, a binary scattering event in homogeneous space is translational 
invariant. As a consequence, one obtains the separation of the center-of-mass motion from 
the relative dynamics, hence a significant simplification of the problem. In optical lattices 
in contrast, we have to use the eigenstates of a lattice, i.e., Bloch-states \q,n) with a quasi- 
momentum q and band-index n. Due to the periodic potential, we find a structured energy 
continuum eiiq, n), which consist of several bands n of increasing width (defining upper and 
lower band edges), as well as forbidden gaps. The density of states also varies accordingly. 

The discrete translational symmetry of the lattice can be exploited for the scattering 
problem by introducing a crystal momentum basis [|23l . which is formed by two-particle free 
Bloch states \qi,ni)i^\q2,n2). They are labeled by the center-of-mass momentum Q = qi + q2 
modulo the crystal momentum fc^, which remains a good quantum number. In this way, the 
coupled two-particle scattering simplifies again and we can introduce a scattering amplitude 
to measure the strength of the binary interaction. 

By introducing the Feshbach coupling, in addition to the normal pairing potential 
between two fermions, we will be able to enhance the strength of binary interactions via 
an external magnetic field [|24l . Being in a lattice, we can consider interband interactions in 
addition to the intraband transitions. This will be important if the strength of the Feshbach 
coupling is comparable to the interband separation. For simplicity, we will consider here 
only very narrow resonances, such that a single band description is sufficient. Extensions to 
multi-band configurations are straight forward in the present formulation, if necessary fl25l . 

In this article, we present the principle of a scattering calculation for two fermions in a 
one-dimensional lattice. In Sec.[2l we will start from the current form of the many -body theory 
that is pursued by many groups. We will simplify this by considering only two fermions and 
a compound bosonic molecule to obtain a two-component two-particle Schrodinger equation 
for the molecule and fermion-pair wave-function. In Sec. [3l we will briefly review the basic 
concept of the scattering phase, the transmission probability and discuss the simplest model 
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for a Feshbach resonance in homogeneous space. This will be generalized to two-particle 
Feshbach scattering in a one-dimensional lattice in Sec HI We present numerical results for 
Feshbach resonance within a single band approximation and demonstrate that it can be tuned 
selectively with the center-of-mass momentum, due to the non-parabolic energy-band. 



2. Reducing many-body physics to two atoms and a molecule 



Currently, much effort is devoted to the many-body description of resonance superfluidity 
and the BEC-BCS crossover |[26l l27l l28ll . Thus, we will briefly introduce the fundamental 
model Hamiltonian ll29l [30l . Then, we will apply it to the situation of only two fermions 
and a bosonic molecule, now, trapped in the same periodic one-dimensional optical lattices. 
The periodic trapping of both molecular and atomic components is beneficial to the overall 
interaction cross section as both components will tend to be localized at the anti-nodes of the 
optical potential, thus be constantly available for scattering. 

In the language of second quantization, we describe the many-body system with 
fermionic fields %Ij^{x), which remove a single fermionic particle from position x in internal 
state a = {t, |}, and molecular bosonic fields (j)i{x), which annihilate a composite bound 
two-particle excitation from the center-of-mass space-point x in internal configuration i. 
These field operators and their adjoints satisfy the usual fermionic anti-commutation rules 



6{xi - X2) 60 



0, 



(1) 



and bosonic commutation rules 



X2 



S{Xi - X2)Si^ 



«2 5 



0, 



(2) 



respectively. We want to consider only bound molecular excitations by choosing a high 
dissociation threshold energy. Effectively this closes this decay channel and allows only 
for collision induced processes. These are the basic ingredients of a Feshbach resonances 
as originally invented by E. Fermi, U. Fano and H. Feshbach [|3T1 [32l [33l . Hence, we can 
present the composite molecular field also with respect to the individual coordinates (xi, X2) 
of the dimer, i. e. 



(xi,X2) = ^0i((xi +X2)/2)(X2 -Xi|i) 



(3) 



The dynamics of the multi-component gas is governed by a total system Hamiltonian 

H = Hl + V (4) 
Hl = jd'x J2 ^i(x)//L(x,p)V^,(x) 



'^={T,i} 



+ / d X0T(xi,X2)ifL(Xi,Pi,X2,P2)0(Xi,X2) 



(5) 
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t>= j <\^X |'t;q(Xi - X2)0^(Xi,X2)0(Xi,X2) 
+ 5f*(xi - X2)0'''(xi, X2)V'^(x2)V'^(xi) + H.C. 

+ ^;p(xi - X2)^}(xi)Vi}(x2)V'^(x2)^/'^(xi)| . (6) 

It consists of the lattice Hamiltonian Hi and the interactions V between atoms and molecules. 
We assume that the free dynamics of the atoms and molecules is determined by their kinetic 
and potential energy in a quasi- Id optical lattice 

/7i(x,p)= -^V^ + UL{x) + U^{y,z), (7) 
zm 

where x = {x, y, z) and p = —ihV are canonically conjugate variables in the position 
representation, m is the atomic mass and Ul{x) is a one-dimensional optical lattice potential. 
Furthermore, we want to assume that the motion in the perpendicular (y, z) direction is 
effectively frozen out by a tight confinement potential U± {y,z). Supposedly, the lattice energy 
is identical for the fermionic atoms of both kinds and twice that for the molecules i. e., 

i?L(Xi, pi, X2, P2) = Hl{Xi, Pi) + i?L(x2, P2). (8) 

The binary interaction potential Vp accounts for the non-resonant interaction of "spin-up" and 
"spin-down" fermions, the coupling strength g converts free fermionic particles into bound 
bosonic molecular excitations and Vg is the molecular potential with at least one bound state. 
One can further simplify matters by considering only even parity binary interactions potentials 
Vg, Vp, and g. Moreover, we have neglected the interactions among the molecules, since we 
will focus on the case of just two fermions. 

The essence of the resonant scattering physics in this many-body Hamiltonian can be 
brought out most clearly, if we consider only two interacting fermionic atoms with field- 
amplitude t/j and a bosonic molecule with amplitude cj), i. e., 

y"d6x|0(Xi,X2,t)0^(Xi,X2) 

+^(Xi,X2,t)l^}(Xi)l^}(x2)}|0). (9) 

According to the Pauli principle, the bosonic and fermionic part of the total wave function 
must be symmetric and anti-symmetric under particle exchange. Thus, if we limit the 
discussion to the fermionic singlet channel, we need to have a symmetric spatial amplitude 
^(xi,X2,i) = V(x2,xi,i), as well as a symmetric molecular wavefunction 0(x2,xi,i) = 

0(Xi,X2,i). 

In this restricted few-particle Fock space with the imposed constraints on the interaction 
channels, finally one finds a two-component Schrodinger equation for the state vector 

X(Xi,X2,i) = ((/)(Xi,X2,i),'0(Xi,X2,i))^ 

i^tX = [Hl(^i, Pi, X2, P2) O 1 + V{xi - X2)] x, 
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In the following, we will also abandon the three-dimensional character of the problem and 
restrict the discussion to a quasi one-dimensional situation, when all dynamics takes place 
along the x-direction. 



3. Scattering in free space 

The previously introduced two-component Hamiltonian contains the essential ingredients of 
Feshbach scattering, but lacks the full translational invariance of homogeneous space, due 
to the presence of the lattice potential. We will therfore briefly review the scattering phase 
and the prerequisites for the appearance of a Feshbach resonance in free space with a contact 
potential for later reference. 



3.1. Scalar potential scattering 

The energy of two particles on a line 

H = Ho{p,) + Hoip2) + V{xi-X2), (11) 

consists of kinetic energy Hq = p'^/2m and the short-range binary potential energy V. 
Obviously, all parts of this Hamiltonian are translational invariant. This symmetry can be 
exploited by introducing a center-of mass coordinate X and a relative coordinate x, as well as 
their conjugate momenta P, p as 

X={xi+X2)/2, X = Xi-X2i (12) 

P=Pi+P2, p={pi-p2)/2. (13) 

The total momentum P is the conserved quantity that is associated with the symmetry 
generating displacement operator (a) = exp [iaP/h]. It shifts the whole system by a 
distance a in the ei2 = ei + 62 direction and commutes with the Hamiltonian 

[H,Vu{a)] = 0. (14) 

Thus, the total two-body wave-function X2{xi,X2) can be expressed in terms of center-of- 
mass and relative coordinates with aproduct ansatz X2{xi, X2) = 9{X)x{x). Here, we assume 
that 9{X) = exp (iQX) is a plane wave in the center-of-mass coordinate with momentum 
Q and the relative wave-function is denoted by xi^)- This reduces the corresponding two- 
particle Schrodinger equation with total energy E to the standard form Q of an effective 
single particle problem 

^ + e-V{x)\xix) = 0, (15) 

where we have also rescaled all dimensional quantities, like length x — > xk^, in terms of the 
wave-number fc^. of an optical lattice photon [see Eq. (|25])1, as well as energies or potentials 
V V/el, in terms of the recoil energy el = h^k\/m. Adopting such units, will facilitate 
the comparison with the lattice case discussed in Sec. HI Then, the relative energy e becomes 
synonymous with the free space dispersion relation 

n2 

e{k) = e = E --^ > 0. (16) 
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As any second order differential equation, Eq. (fTSl ) admits two linearly independent 
scattering solutions Xk^\^) for the same energy e{k). By comparing them to non-interacting 
plane waves, one can introduce reflection and transmission amplitudes R(k) and T{k), i. e., 

lim x^^\x,k) = e'^"" + Re-'^'' 

a;— »— oo 

= e'^" cos{kx - 6e) + ie'^° sm{kx - 60), (17) 

X— >oo 

= e'^" cos{kx + 6e) + ie'^° sm{kx + 60) (18) 

as well as even and odd scattering phase shifts 5e{k) and 5o{k). The transmission and 
reflection probabilities are given by 

|T|2 = cos2(5, - 60), \R\^ = sm\6, - (5„), (19) 

and unitarity demands current conservation, which is mathematically paraphrased as |Tp + 

= 1. 
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Figure 1. Even scattering phase Se{k) (solid line) and transmission probability |Tp (dashed 
line) of a Feshbach resonance as a function of energy e — k"^ with a resonance at e^es — 6. 
Parameters: m = 10, w = —4 and g — 1. 




3.2. Two-component Feshbach resonances with a contact potential 

The quintessential mechanism for a Feshbach resonance occurs, when one couples two 
internal states x = iXgi Xp)~^ to the external motion of the one-dimensional, two-particle 
Hamiltonian of Eq. (fTTI) . The closed q-channel needs an attractive delta potential to support a 
bound state and a localized coupling matrix element is required to provide an interaction with 
the open p-channel [|34l . A simple, analytically solvable model for this is 

H = Ho{pi,P2) ® 1 + V{xi - X2), 
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V{x) 



u + v6{x) gS{x) 
g5{x) 



(20) 



The parameters of the potential matrix V are the threshold energy m > 0, the strength of 
the closed channel potential < and the channel coupling parameter g. As in the scalar 
case, all contributions to the Hamiltonian are translationally invariant and a product ansatz 
for the wave-function X2(^i5 ^2) = exp {iQX)x{x) leads to the Schrodinger equation for the 
relative wave-function 



0. 



(21) 



Here, we have again introduced rescaled potentials, energies and length mentioned in the 
context of Eq. (fT5l) . 

For relative scattering energies < e{k) = k"^ < u, smaller than the threshold energy 
u, the closed channel asymptotic wave-function Xq{^) vanishes exponentially, while the open 
channel wave-function propagates outward freely. Due to the even parity of the contact 
potential, we can also choose wave-functions of definite parity, i. e.. 



lim Xq^^ (^) k) e 



-q\x\ 



>itoo 



lim Xp\x,k) 



q{k) 

cos{kx ± 6e). 



Vu — k'^, 



(22) 
(23) 



Odd wavefunctions do not accrue any phase shift, as they vanish at the point of interaction. 




Figure 2. Even scattering phase Se{k{Q, E)) versus energy E and momentum Q. The 
Feshbach resonances are located on a free space parabola as a function of Q. Parameters 
as in Eq. (|32] |. 



By integrating the coupled Schrodinger equation Eq. (|2TI) over an infinitesimal strip 
around the discontinuity, one obtains the necessary boundary conditions to compute the 
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scattering phase of the even wavefunction. After some minor algebra, one obtains the phase 
shifts 



This phase shift and the corresponding transmission probability |Tp of Eq. (fT9l) are depicted 
in Fig. [TJ One clearly observes the Feshbach resonance with a vr-phase jump at er^es = 
fc^gg = u — The position of the resonance energy is usually controlled by changing 

the dissociation threshold energy u via magnetic fields or via the depth of the bound state 
V, which is an property of the considered atomic element and, thus, harder to modify. The 
width of the resonance is primarily determined by the coupling strength g ll35l . Due to the 
vanishing odd phase-shift, the total reflection T = also occurs at the same energy, when the 
even phase shift reaches SeiKes) = 7r/2. It is also important to note that the phase can be well 
approximated by a Breit-Wigner curve in the vicinity of the resonance. 

In Fig.O we present the same scattering phase 5e{k{E, Q)) of Eq. ((24)) . now as a function 
of the two-particle energy E and the center-of-mass momentum Q, according to Eq. (fT6l) . This 

representation emphasizes the parabolic form of the resonance energy ii^res = > 0, 

as a function of Q. The Feshbach resonance is always present for e^es > and for all Q, or it 
disappears if Sres < 0. As we will show in the following sections, this fundamental behavior 
can be changed by considering Feshbach scattering in a lattice where center-of-mass and 
relative motion are coupled. Results can be seen in the final Figs. [TTI and [T2| 



4. Scattering on a one-dimensional lattice 

4.1. Block states and periodic boundary conditions 

Studying physics in periodic structures immediately leads to the consideration of Bloch 
states ll36l l2ll. which reflect the discrete translation symmetry V{la) = exp [ilap/h] of the 
Hamiltonian [HL,V{la)] = with / G IN . We have already introduced a quasi-one 
dimensional lattice Hamiltonian in Eq. ([7]) and will disregard from now on the transverse 
degrees motion, i. e., 

Hl{x,p) = ^ + Ul{x), Ul{x) = Ulcos^ {kLx/2). (25) 

The lattice potential f/z, (x + a) = f/^ (x) is characterized by a lattice constant a, which in turn 
defines the crystal momentum ki = 27i/a. Thus, the eigenstates of Schrodinger equations 

HL\q,n)=el{q)\q,n), (26) 

can be classified as Bloch states \q, n) with quasi momentum — vr/a < q < vr/a, and energy 
bands s'^{q), labeled by band-index n. According to the Bloch theorem 

(x|g,n) = e^^XW' <(a; + a) = (27) 
Such an energy eigenfunction can be decomposed into a plane wave phase factor and a 
lattice periodic function vl^. Subjected to translation to the next lattice site, the wave-function 
acquires a complex phase 

{x + a\q,n) = e''^''{x\q,n). (28) 
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The periodic continuation in momentum space has its subtleties and attention needs to be paid 
to the vanishing of the wave-function at the center or edge of the first Brillouin zone [|36ll . 
However, the generic case is simply given by {x\q + ki^n) = {x\q,n). 

0.8 
_ 0.6 
0.4 
0.2 


-1.5 -1 -0.5 0.5 1 1.5 

X 





Figure 3. Magnitude (a) and phase 'y{x) (b) of the lowest band Bloch function 

{x \ q,n = 0) — \u]^{x)\exp[ij{x)], versus dimensionless position —tt/2 < x < tt/2, 
for a few momenta q ~ (—1, —0.5, 0, 0.5, 0.992). In this shallow potential the Bloch waves 
are almost plane and deviations are only seen at the band edge q « ±1. Parameters as in 
Eq. (O. 



In practice, it is usually necessary to work with a discrete subset of Bloch states, which 
are found by considering a periodically continued, finite lattice with an even number of wells 
= 2M. The Born-von-Karman periodic boundary conditions, then require that 

{x + L\q,n) = e*'"^ {x \q,n) = {x\ q, n) . (29) 

where the total length L = Na. This can only be true, if there are exactly N distinct momenta 
qi in each band 

TT 27r/ / TT , 

--<qi = -r = kL^<-. -^<l<^- (30) 
a L JM a 2 2 

We choose the following orthogonalization for the Bloch-states on a finite lattice 

{qi,ni\q2,n2) = Ndq^^q^dn^^ni, / dx\{x\q,n)\'^ = 1. (31) 
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The general behavior of the Bloch states is shown in Fig. [3l There we have selected a 
few Bloch eigenfunctions for some momenta q within the Brillouin zone, which exemplify the 
modified plane wave behavior. A typical energy-band structure in coordinate and momentum 
space is presented in Fig. |4] and Fig. [51 respectively. In here and all of the following 
calculations, we use dimensionless parameters for a shallow lattice that only supports one 
band below the barrier i. e., 

Ul = -(6/5)^ a = TT, = 2^ (32) 





Figure 5. Energy bands e° (9) (solid) and e\{q) (dashed) versus quasi-momentum q in the 
first Brillouin zone —l<q<l. 
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4.2. Scalar two-particle scattering in lattices 

We will now consider the scattering of two particles in the presence of a lattice. Thus, we 
have a two-particle lattice Hamiltonian and a binary interaction V 

H = Hl{xi,Pi) + Hl{x2,P2) + V{xi - X2). (33) 

By forming two-particle basis states out of single particle Bloch waves one obtains two- 
dimensional Bloch states 

|Q,g,ni,n2) = |gi(<5, g), ^1) ® |g2(<5, g), ^2) (34) 

Q = qi + q2, q = (qi - g2)/2 (35) 

qi = Q/2 + q, q2 = Q/2 - q. (36) 




Figure 6. Reciprocal momentum space for two particles in a one-dimensional lattice 
consisting of only = 4 wells. The Cartesian axes represent the individual momenta in 
the reduced zone scheme —1 < qi,q2 < 1- The 45° rotated grid lines correspond to the 
contours of center-of-mass — 1 < Q < 1 and relative momentum — 1 < (? < 1. 

The corresponding reciprocal momentum space for two particles in a one-dimensional 
lattice is depicted in Fig. |6l One can see the individual single particle momenta qi and q2, 
as well as the center-of-mass momentum Q and relative momentum q. This set of symmetry 
adjusted basis states can be used to represent a general quantum state with quasi-momentum 
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Q as 

\x.Q)= E (37) 

i},ni,n2 

^^i2Mlx,Q) = e*'"«|x,Q). (38) 
From the discrete translation symmetry of the system, one finds the selection rule 

= {x,Q\[H,V,,{la)]\x,Q') = {x,Q\H\x,Q'){e''^^'^'-'^^ - 1), (39) 
which implies Q' = Qmod/c^. In the following, we will consider only transitions within 

the lowest energy band i.e., rii = n2 = and simplify the notation to Xq''^\(1) = Xq(Q')' 
consequently. 




Figure 7. Energy contours El {Q, q) for the non-interacting two-particle system in momentum 
space. The superimposed Q =const. line (solid black) has two intersections with another 
E —const, contour (dark gray), which are at q (circled) and —q (boxed). 

Within these assumptions, one can represent the Schrodinger equation to Eq. (l33l) as 

0=J2{i'Liqi)+eliq2)-E]5,,, + VQiq,q')}xQiq'), (40) 
<?' 

where qi = qi{Q, q) and q2 = q2{Q, q)- If we let our perturbation tend to zero, then we are 
left with the two-dimensional energy surface 

^L(g,g) = 4fe) + 4fe). (41) 

The contour lines of which are depicted in Fig. |71 The marked intersections show that there 
are two linearly independent two-particle quantum states labeled by (gi, 52) and (g2, qi) that 
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have the same total energy E and center-of-mass momentum Q. It it interesting to note that 
in general there are no straight contour- lines atQ = ±1 or g = ±1, as would be the case in 
the tight binding limit \Ul\ ^ 1- The energy range that is covered by the two-particle energy 
can be seen in Fig.[8l It shows the projections along the Q and q momentum lines. 
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Figure 8. Projections of the two-dimensional energy surface El{Q, q) versus momentum Q 
(a) and q (b). 





Figure 9. Interaction matrix element VQ{q,q') in reciprocal momentum space for Q 
other parameters as in Eq. ( [32] l. 



Now, if we turn to the evaluation of the matrix element Vgi^q, q'), then it is best expressed 
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in center-of-mass coordinates. Moreover, the algebra simplifies considerable, if we use a zero- 
range contact potential for the binary potential V{y) = v5{y) and obtain 

V{Q,q,(^) = vvQ{q,q'), (42) 

'VQ{q,q') = -j;^[u*Q/2+q{x)u*Q/2^g{x)UQ/2+q'{x)UQ/2-q'{x)y (43) 

The shape of this potential matrix can be calculated quite accurately with localized Wannier 
functions in the Gaussian approximations. In general, this agrees well with the numerical 
results shown in Fig. |9l In this picture, we have chosen a value of Q = for the center-of- 
mass momentum. For other values Q 0, one obtains a modest variation of the shape of the 
matrix element at the momentum edges, but it remains predominantly constant inside. 

4.3. Two component Feshbach scattering in lattices 

Having established the basic notions and concepts in the previous sections, we can now turn 
to the Feshbach resonance scattering phenomenon in a lattice. The basic Hamiltonians have 
been already introduced in Eqs. (fTOl) . (|20l) and (l33l) . Explicitly, we have the motion in the 
lattice Hamiltonian Hi, as well as the Feshbach potential V 

H = Hl{Xi,Pi,X2,P2) ® 1 + V{xi - X2) 

Here we have deliberately disregarded the pairing potential Vp = 0, as it only provides a 
modification to the Feshbach phenomenon. The state of the coupled two-component system 
in the lowest Bloch band n = is now defined as 

ix,Q)^i:(::i:i)i«-). 

where we have characterized the state with the momentum Q and labeled the open and closed 
channels by the molecular </> and two-fermion wave-function tjj as in Eq. (flOl) . 

In order to obtain the Bloch representation of the Schrodinger equation to Eq. (|44)) 



E\x,Q) = H\x,Q), (46) 
we project it on the two-particle Bloch eigenstate and get 

q' 

o))(':i''l)- 

We have numerically diagonalized this Schrodinger equation for the lattice parameters 
of Eq. (|32|) and the Feshbach parameters u = 1.5, v = —7, g = 0.03. Due to the parity 
of the Hamiltonian, one can sort the resulting 2N eigenstates into even and odd states and 
assign them a positive or negative quantum index — < u < N.By turning off the coupling 
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between the manifolds i.e., g = 0, one can distinguish easily the eigenvalues Ey{Q,g = 0) 
that belong to molecular (q-channel) or the open fermionic spectrum (p-channel). We have 
depicted a few representative values in Fig.fTOk). It can be seen clearly that there is a single 
bound state embedded in the allowed energy band of the open p-channel. In Fig. [TOb). we 
present the interacting eigenvalues Ey{Q, g = 0.03), which are now formed by admixtures of 
both manifolds and can only be distinguished by the parity of the state. 
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Figure 10. Discrete eigenenergies E^{Q,g) versus quantum number —N < v < N for 
momentum Q — 0.41. Negative quantum numbers label odd states, while positive numbers 
count even states, a) Energies of uncoupled system E^{Q, g ~ 0). Note the location of the 
bound state (*) of the molecular manifold (0, □) is inside the fermionic spectrum marked by 
(x, o). b) All interacting eigenenergies E^{Q, g = 0.03) versus quantum index i'. Odd states 
(solid) and even states (dashed dotted). The bound state is now an admixture of molecular and 
fermionic manifolds. 



In order to describe the scattering resonance, we use the standing wave version of the 
Lippman-Schwinger equation Q 

\X^'\ Q, g) = ep ® \Q, q) + G^I\e)V\x^'\ Q, q), (48) 

where Gp ® \Q,q) represents a non-interacting Bloch-wave in the asymptotically open p- 
channel Bp = (0, 1)^. For the lattice Hamiltonian H^, one obtains a standing wave Greens 
function G^l^ from the usual retarded and advanced Greens functions G^^^ by 

Gt\E) =\{G^^\e) + g[^\e)), (49) 

G'^^\e) = lim . (50) 

L V ; e^o+ E±ie- Hl 

Finally, the effect of scattering is measured by the off-shell Heitler- or K-matrix element 

KQ{q',q) = {Q,q'\eJV\x^^\Q,q). (51) 
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We assign an even on-energy shell scattering phase shift 

tan [6,{Q, E)] = Kgiq, q'), E = E{Q, q) = E{Q, q'). (52) 

We have evaluated this phase shift for two differently strong bound molecular states and 
varied the center-of-mass momentum of the two-particle collision. In the first Fig.[TTl one can 
see the Feshbach resonance cutting into the energy band and disappearing once it is outside. 
This is modified in Fig. [I2l which has a lesser bound resonance, and therfore does not leave 
at the lower band edge. 

This has to be compared with the Feshbach scattering in homogeneous space according 
to Eq. (l24l) and has been presented already in Fig. [21 There, the resonance energy E = 
£res + <5^/4 > 0, is on a simple parabola as a function of Q. If we tune the Feshbach 
resonance energy Sres < 0, it either disappears completely or is always present. 




Figure 11. Even scattering phase-shift 5e{E, Q) of a two-particle Bloch-wave versus energy 
E and momentum Q. Once can see clearly Feshbach resonances symmetrically located around 
the Q-axis, which disappears when the bound state is outside the band. Parameters; u = 1.5, 
w = -7,3 = 0.03. 



5. Experimental considerations 

Starting from a dilute Fermi gas, atoms can be loaded adiabatically into the lowest energy band 
of the lattice. For quasi-momenta not too close to the band edge, this requires that the Fermi 
energy remains clearly below the atomic recoil energy ^r2_^corresponding to In ■ 8.4 kHz 
and l-K ■ 74 kHz for the and ^Li D2-lines respectively |j. To ensure that the Feshbach 
resonance provides coupling only within a single energy band, its width has to be sufficiently 

I Here hiOr — ft^fc^/2m = 6^/2 represents half the energy unit defined Eq. (fTsT l. 
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narrow, i. e., below a value of about huj^/ /S.ji, where A/i denotes the difference of the atomic 
and the molecular states magnetic moments respectively. For typical parameters, we arrive 
at a required width of the Feshbach resonance below 10 to 100 mG. We are aware, that the 
investigation of such weak Feshbach resonances requires a high magnetic field stability. On 
the other hand, three-body losses, which for free space experiments are a mayor second 
experimental constraint in the study of narrow Feshbach resonances, are expected to be 
reduced in the presence of the lattice potential due to atom localization in the micropotentials. 

For an experimental verification of the relative atomic momentum dependence of the 
Feshbach resonance, as indicated in Figs. [TT] and [121 one could study such scattering 
resonances for a variable filling of the lowest energy band in the lattice. In this way, due 
to the Pauli exclusion principle, for an increased band filling different regions of the relative 
atomic quasi-momentum are subsequently filled up. The variation of the lineshape of the 
Feshbach resonance is determined by the coUisional properties of atoms in the lattice. For 
example, if we take the situation of Fig. [TTl for small atom filling the scattering resonance 
would disappear. The Feshbach resonance here could only be observed for a large atom 
filling. 



1^ 




Figure 12. Even scattering phase-shift 5e{E, Q) of a two-particle Bloch-wave versus energy 
E and momentum Q. With a less deeply bound state the resonances remain above the lower 
band edge. Parameter; u — 1.5, v — —5, g — 0.03. 

6. Conclusions and Outlook 

We have considered Feshbach scattering of two Fermions in a one-dimensional optical 
lattice. Due to the reduction of translation invariance, the relative and center-of-mass motion 
are coupled and we present the scattering theory in the crystal momentum basis. Within 
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a single band approximation, we have calculated numerically Feshbach resonances and 
demonstrate that the position of the Feshbach resonance depends selectively on the center-of- 
mass momentum, due to the non-parabolic shape of the energy band. The simple resonance 
structure suggests that a semi-analytical calculation based on an effective-mass approximation 
should be applicable and is currently pursued. 

In the present article, we have only considered Feshbach resonances that are much 
weaker than the interband separation. Thus, we could limit the discussion to a single energy 
band. Wider Feshbach resonances will involve the coupling of several bands. Going beyond 
the single band approximations will provide interesting new physics as already advanced in 
Ref. [25] for the case of potential scattering. This will be examined, based on the present 
formulation, in a forthcoming publication. 
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